Efficient Linear Optics Quantum Computation 
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We investigate the computational power of passive and active linear optical elements and photo-detectors. 
We show that single photon sources, passive linear optics and photo-detectors are sufficient for implementing 
reliable quantum algorithms. Feedback from the detectors to the optical elements is required for this implemen- 
tation. Without feedback, non-deterministic quantum computation is possible. A single photon source sufficient 
for quantum computation can be built with an active linear optical element (squeezer) and a photo-detector. The 
overheads associated with using only linear optics appear to be sufficiently low to make quantum computation 
based on our proposal a viable alternative. 



I. INTRODUCTION 



Quantum computers promise to greatly increase the effi- 
ciency of solving problems such as factoring large integers [|l]j 
and combinatorial optimization [0]. One of the greatest chal- 
lenges now is to implement the basic quantum-computational 
elements in a physical system and demonstrate that they can 
be reliably and scalably controlled. One of the earliest pro- 
posals [J3J for implementing quantum computation is based on 
encoding each qubit in two optical modes together containing 
exactly one photon. The main problem with this proposal is 
that it is extremely difficult to non-linearly couple two optical 
modes containing few photons. Here we consider the question 
of what can be accomplished in principle using various com- 
binations of only the simplest optical elements: passive linear 
optics, photo-detectors, squeezers (an active linear optical ele- 
ment) and non-deterministic single photon sources. We show 
the surprising result that quantum computation is possible in 
principle using the first three of these. Only weak squeezing 
is required, and this only for single photon generation. Al- 
ternatively, any non-deterministic single photon source can be 
used. The optical elements must be controllable based on out- 
put from the photo-detectors. An alternative method relying 
on phase space based encodings, homodyne detection and an 
optical non-linearity in state preparation was independently 
presented by Gottesman and Preskill [|J. This establishes 
the principle of linear optics quantum computation (LOQC). 
Furthermore, the basic elements are testable in today's lab- 
oratories. Our non-deterministic optical gates may find im- 
mediate application in quantum communication, experiments 
with entanglement and optical state preparation. Spin-offs 
include near-deterministic quantum teleportation and a par- 
ity measurement for photon based qubits. Scalable quan- 
tum computation using our methods requires highly efficient 
photo-detectors, very low loss short term photon storage and 
long state preparation times to achieve the minimum accu- 
racy required for reliable quantum computation. However, 
the theoretical overheads are sufficiently small compared to 



the implementation-independent requirements for scalability 
to suggest that LOQC is a viable proposal for implementing 
quantum computers. 

We begin by briefly introducing the basic notions of 
bosonic modes and linear optics; give the representation of 
qubits using bosonic modes; describe the basic techniques re- 
quired for implementing non-deterministic and and then de- 
terministic quantum computation. We then consider scalabil- 
ity issues and argue that using quantum codes, the complexity 
of LOQC can be significantly reduced. We conclude with a 
discussion of the practical issues involved in using our meth- 
ods. We assume some familiarity with quantum optics and 
quantum computation. An introduction to quantum optics can 
be found in [Q]. An overview of quantum computation and 
implementation issues is in | 



II. BOSONIC MODES 

To avoid issues surrounding the polarization states of pho- 
tons we cast our work in terms of non-interacting spin-less 
bosons. Thus, the physical system of interest consists of a 
number of bosonic modes. A bosonic mode (or quantum har- 
monic oscillator) is a quantum system whose state space is 
spanned by the number states |0), |1), |2), . . .. The modes are 
labeled and, with explicit labels, a number state of mode / 
is denoted by \k\. The vacuum state, which satisfies that all 
modes are in state |0), is denoted by |0). The observables 
for mode I can be constructed from the annihilation opera- 
tor sS l \ which satisfies a^|fc^ = \f~k\k — l\ for k > 1 and 

a^|0) = 0. and n' 1 ' := a^ a^ are the creation and 
the number operator for this mode, respectively. 



III. OPTICAL ELEMENTS 

The most readily implementable processes are those given 
by passive linear optical elements. These are elements whose 
effect on the state of the mode are given by the following two 
Hamiltonians: 



knill@lanl.gov 



I aflamme @ lanl. go\ 



riilburn@physics.uq.edu.au 



n 

h (lr. 



(/) 



a (0t a (0 



(phase shifter) 
a^aM 1 + ad'a^) (beam splitter). 
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Evolutions implementable by passive linear optics preserve 
the total number of bosons in the modes. It is therefore con- 
venient to describe them by their effect on the creation opera- 
tors. Specifically, if U is the unitary operator associated with 

the evolution, then U takes the state a^ T |0) to 



L/a w V|0) 



(2) 
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using the fact that U'\0) = |0). The matrix defined by the 
coefficients U m i must be unitary, and furthermore, for all uni- 
tary U m i, there is a sequence of phase shifter and beam split- 
ter evolutions that implements the corresponding operation up 
to a global phase For a named optical element X, let 
U(X) be the unitary matrix associated with X according to 
the above rules. The unitary matrices associated with phase 



shifters P s ^ l> and beam splitters B g {lm> are: 



(Zm) 



U{P e {1) ) = e 



cos(0) — sin(0) 
sin(0) cos(0) 



General linear optical elements have Hamiltonians that con- 
sist of terms at most quadratic in the annihilation and creation 
operators. We will make use of the squeezer, an active linear 
element whose Hamiltonian is 
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A mode can be prepared in its vacuum state |0). A non- 
deterministic single boson source prepares a given mode in 
the state |1) with some probability of success. It is assumed 
that the information on whether the source was successful is 
available classically. 

Measurement is accomplished by means of a particle de- 
tector, which destructively determines whether one or more 
bosons were present in a mode. We assume that measurement 
results can be used to control other optical elements. A more 
powerful detector is a destructive particle counter, which re- 
turns the number of bosons present in a mode. An approx- 
imate particle counter that suffices for our purposes can be 
designed by using N particle detectors operating on N modes 
ll, . . . , In. To measure mode I, use passive linear optics to 



transform a' 1 ' 



l 



,(I»)t 



Each of the N modes 



is then measured using the corresponding particle detector. If 
mode I has k bosons, then the probability that some output 
mode has two or more bosons is at most 1 



(AT) 



N k 



< fc(fc-l) 
2N ■ 



where we define (N)k as the fc'th falling factorial of N, given 
by N(N - 1) . . . (N — k + 1). Provided that the maximum 
number of photons is not too large, the number of bosons is, 
with high probability, the number of particle detectors that de- 
tect a boson. The bound on the error probability is required for 
estimating the resources needed to achieve sufficiently high 
accuracy for reliable quantum computation. 



IV. REPRESENTING QUBITS 

A qubit is a quantum system with two distinguished basis 
states, |0) and |1) . Thus, the state space of a qubit general- 
izes that of a classical bit by allowing superpositions of the 
classical states. To represent a qubit using bosonic modes 
requires mapping the basis states unitarily into two states of 
one or more modes. We choose the traditional encoding us- 
ing two modes and one boson, where |0) — ► lO^Jl), and 
|1) — > |1)J0^, and call this a "bosonic qubit". For exam- 
ple, the two modes can be two orthogonal polarization states 
of an optical mode, so that the polarization state of a photon 
encodes the state of a qubit. With this encoding, the one qubit 
(U(2)) rotations are easy to implement using passive linear 
optics. In the case of encoding in polarization states, polariza- 
tion rotators can be used in place of beam splitters. 

To complete the implementation of a standard quantum 
computer, it is necessary be able to prepare the state of a qubit 
in |0) , to measure a qubit in the standard basis \0) q , |1) , and, 
in addition to the rotations above, to implement a controlled- 
not (c-not) and/or a controlled-sign (c-sign) gate [^[. The 
state preparation is accomplished here by using a single bo- 
son source to prepare mode a in |l) a and mode b in \0\. To 
measure a bosonic qubit it suffices to use a photo-detector on 
the first mode. Unfortunately, linear optics is insufficient for 
implementing the two qubit gates c-not and c-sign, which are 
defined in terms of unitary matrices acting on the four dimen- 
sional two-qubit state space by 



c-not 



c-sign 
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where the ordering |00), 1 01) , 1 10) , 1 1 1) is used for the basis. 
Thus c-not flips the second qubit conditional on the first, and 
c-sign flips the sign of the second qubit conditional on the first. 
A direct bosonic implementation using a suitable Hamiltonian 
requires strong optical non-linearities. For example the c-sign 
gate could be implemented with a non-linear phase shift on the 
second halves of two bosonic qubits. The task of this work is 
to show how these operations can be induced using only single 
particle sources and particle detectors. 

We conclude this section by introducing one more useful 
gate, the Hadamard transform on one qubit, defined by the 
matrix 



H 



1 

71 



i i 

i -i 



(7) 



H can be implemented by a balanced beam splitter on bosonic 
qubits (to within a global phase). Note that c-not can be ob- 
tained from c-sign by conjugating the second qubit with the 
Hadamard transform (which is its own inverse). 
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V. NON-DETERMINISTIC QUANTUM COMPUTATION 



VI. NEAR-DETERMINISTIC QUANTUM COMPUTATION 



A non-deterministic operation is one that succeeds only 
with probability < 1, and whether it succeeded or not is 
known after it has been applied. It turns out that it is much 
easier to implement a non-deterministic non-linear phase shift 
rather than a deterministic one. These non-deterministic 
gates will also play a crucial role in preparing the states re- 
quired for deterministic non-linear phase shifts. At the same 
time, having such a gate implies the ability to perform non- 
deterministic quantum computation (conditional on an accept- 
able measurement outcome at the end), general state prepara- 
tion, and entanglement sharing through lossy systems. The 
basic nondeterministic gate we construct performs the trans- 
formation 

NSi : oo|0) +ai|l) + a 2 \2) -> a |0) + a x |l) - a 2 \2) 

with probability 1/4. (8) 

In order to implement the conditional sign flip c-sign on two 
bosonic qubits encoded in modes 1,2 and 3,4, respectively 
do the following: Apply the balanced beam splitter B jr ^ 2 to 
modes 1 and 3. In the case where the two qubits are both in 
state |1^, modes 1 and 3 are in the state |11\ 3 , which trans- 
forms to -^(120^3 + 102)^). In none of the other cases do two 
bosons appear in the same mode. Now apply NS\ to mode 1 
and then to mode 3 and undo the balanced beam splitter. This 
has the desired effect with probability 1/16. 

The procedure for applying NS\ to mode 1 begins with 
preparing the initial state |10), 3 in ancilla modes 2 and 3. (The 
labels are chosen for convenience and are arbitrary.) Next a 
sequence of beam splitters is used to implement the symmet- 
ric unitary transformation 

/1-2V2 2-V* (3/2 1 /2_ 2 ) 1 /2 

U= 2- 1 / 4 1/2 1/2-1/2 1 / 2 

V (3/2 1 / 2 - 2)V2 1/2 -1/2V2 2 1 /2_i/ 2 

(9) 

Finally, modes 2 and 3 are measured, and the outcome is ac- 
cepted only when the state is found to be |10^ 3 , that is, the 
state is the same as the initial ancilla state. To check that this 
works, observe that due to particle conservation, the condi- 
tional state transformation is of the form 

a |0)+ai|l)+a 2 |2) -> A a |0) + Aiai|l)+A 2 a 2 |2). (10) 

The Afc are given by 

Ao = U22 

= 1/2 (11) 

Ai = U11U22 + U12U21 
= (l-2 1 / 2 )/2 + 2" 1 / 2 

= 1/2 (12) 
A 2 = U 11 (U 11 U 22 + 2U l2 U2i) 

= UuiUuUn + UuUiJ + UuUuUii 
= (l-2 1 / 2 )/2 + (l-2 1 / 2 )2- 1 / 2 
= -1/2 (13) 
A direct calculation shows that U's columns are orthonormal. 



As noted earlier, it is in principle enough to implement a 
deterministic c-sign gate. Due to the results from fault tol- 
erant quantum computing, it is acce ptabl e for the gate to fail 
with a small probability (see Sect. VIII). The basic idea in 
implementing a near-deterministic c-sign gate with low error 
is to use ideas from quantum teleportation || [l0|]. The trick 
is to prepare an appropriate entangled state suitable for tele- 
portation with the desired gate already applied, before using 
it for the teleportation protocol. The problem then becomes 
that of preparing the entangled state (which can be done non- 
deterministically) and implementing the requisite measure- 
ment in the protocol. 



A. Quantum teleportation 

A basic quantum teleportation protocol for transferring the 
state ao\0\ + cei\l\ to mode 3 first adjoins the "entangled" 
ancilla state (\01\ 3 + \10\ 3 )/y/2 to mode 1. Note that in this 
case, the ancilla state is easily generated from 1 10) 23 by means 
of a beam splitter. The second step is to measure modes 1 and 
2 in the basis (|01) 12 ± \10\ 2 )/V2, (|00) 12 ± \ll\ 2 )/V2 (the 
"Bell basis"). We decompose the measurement into two steps. 
The first step determines the parity p of the number of bosons 
in modes 1 and 2 ("parity measurement"). The second deter- 
mines the sign s in the superposition. Consider the case where 
p = 1. Then if s = '+', the state of mode 3 is «o|0) 3 + a i|l^- 
If s = the state is aolO^ — ai|l)j, which can be restored to 
the starting state by using a phase shifter. For p — 0, the situ- 
ation is similar except that 1 0^. and \1\ are flipped (and cannot 
be easily un-flipped using linear optics). The key property of 
quantum teleportation is that the input state appears in mode 3 
up to a simple transformation without having interacted with 
mode 3 after the preparation of the initial ancilla state. 

Consider the parity measurement. Applying a balanced 
beam splitter to modes 1 and 2 and then measuring the num- 
ber of photons in the two modes successfully determines the 
parity, and if it is odd, the sign. As a result, this measurement 
method can be used to perform the teleportation with success 
probability 1/2. We refer to this partial Bell-state measure- 
ment as BMj . 

Two teleportation steps using BMj can be used to imple- 
ment c-sign with success probability 1/4. To see how to do 
this, observe that to implement c-sign on two bosonic qubits in 
modes 1, 2 and 3, 4 respectively, we could try to first teleport 
the second modes of each qubit to two new modes (labelled 6 
and 8) and then apply c-sign to the new modes. The full tele- 
portation procedure requires a correction step that consists of 
applying a Pauli operator (nothing, sign-flip, bit-flip, or their 
product) to each new mode. Since c-sign is in the normalizer 
of the group generated by the Pauli operators, the correction 
step can be deferred until after applying c-sign. Of course, 
now there is nothing preventing us from applying c-sign to the 
prepared entangled state before performing the measurement. 
Thus the implementation of c-sign is now reduced to the prob- 
lem of preparing a modified four mode entangled state 1 64^78 , 
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then implementing the Bell measurement and correcting the 
resulting state in the new modes. The modified entangled state 
is given by 

%) = (|1010) + |0110) + |1001) - |0101))/2. (14) 

| b' 4 ) can be generated with linear optics by following the se- 
quence 

|01)|01) (|01) + |10))(|01) + |10»/2 (with two B„ /2 ) 
\b' A ) (with^i). 

(15) 

The success probability of this procedure is 1/16, which 
means that the expected effort for obtaining one instance of 
| b' 4 ) requires 16 trials of the sequence above. Once such an 
instance is obtained, we can use two BMj measurements fol- 
lowed by phase corrections to implement c-sign on two input 
bosonic qubits with success probability 1/4. To increase the 
success probability further we can either change the telepor- 
tation procedure or improve on the parity measurement. We 
next show how to teleport nearly deterministically and how to 
use this for a near-deterministic c-sign. 

B. Near-deterministic quantum teleportation with beam 
splitters 

An initial entanglement that results in successful teleporta- 
tion with probability of success l/(n + 1) is given by 

n 

\tn) =^|1) J |0)"^|0)^|1)"^ ; (16) 

where we omitted normalization constants. The notation la)-? 
means \a)\a) . . ., j times. Note that if we label the modes 
1 . . . 2n (left to right), the state exists in the space of n bosonic 
qubits, where the fc'th qubit is encoded in modes k and n + k. 
Also, \ti) is the same as the entanglement used in the basic 
teleportation protocol of the previous section. 

The teleportation protocol using state \t n ) teleports mode 
(say) in a superposition of |0) and |1) to one of the last 
n modes of \t n ) by a measurement BM n involving only the 
mode to be teleported and the first n modes of \t n ). The mea- 
surement BM n is implemented using an n + 1 point Fourier 
transform. Let F n be the unitary matrix defined by 

{F n ) u = cj h /V^+T, (17) 

where uj = e l2 ^/( n + 1 ) and k, I £ 0...n. F n is unitary 
and therefore implementable with passive linear optics. Us- 
ing the parallel fast Fourier transform (see page 795 of [pi]]), 
it can be implemented with 0(n\og(n)) elements and depth 
0(log(n)), for n a power of 2. Alternatively, a multiport in- 
terferometer can be used [jj~2|] . Denote the optical array for im- 
plementing F„ by F n . To perform BM n apply F n to modes 
. . . n and measure the number of bosons in each of these 
modes. Suppose we detect k bosons altogether. We claim that 
if < k < n + 1, then the teleported state appears in mode 



n + k and only needs to be corrected by applying a phase shift. 
The modes 2n — I are in state 1 for < I < (n — k) and can 
be reused in future preparations requiring single bosons. The 
modes are in state for n — k < I < n. If k = we learn 
that the input state is measured and projected to |0), and if 
k = n+ 1, it is projected to The probability of these two 
events is l/(n + 1), regardless of the input. We will make use 
of the fact that failure is detected and corresponds to measure- 
ments in the basis |0), |1) with the outcome known. Note that 
both the necessary correction and which mode we teleported 
to are unknown until after the measurement. 

To prove the claim, observe that applying P2-nk/{n+i) to 
mode k for < k < n after applying F n is equivalent to shift- 
ing modes . . . n circularly right before applying F n . This 
means that the states \xi) = F„|l) fc |0)" _fe+1 and \x 2 ) = 
F n \0)\l) k \0) n ~ k differ only by phases in the number basis. 
Thus the measurement BM n cannot distinguish between the 
two states. If the measurement detects rj bosons in mode j 
r j = fc in this case), then the relative phase of the second 
state is given by Ylj ^> r ^ ■ Because of the way \t n ) is entan- 
gled with the last n modes, the measurement outcome trans- 
fers a superposition of \x\) and \x 2 ) in the entanglement to a 
superposition in mode 2n — k + 1. The remaining properties 
of the claim are immediate. 

The teleportation trick to improve the probability of success 
of c-sign can be used with \t n ) also. The necessary modifi- 
cation of \t n \ 2n \t n \ n+1 4 „ to teleport two modes with an 
effective application of c-sign is accomplished by applying 
c-sign to each pair of modes (n + k,3n + I) with k and I in 
1 . . . n. This works up to qubit sign flips because the output 
modes that do not receive the teleported state are in known 
configurations after the procedure. Explicitly, the state that 
needs to be prepared is given by 



K) = (-i) (n ^ )( "~ J) |i) J |o)"- j |o) j '|i)"^ 

i,j=0 

xliyior^lO/'ll)™- 1 . (18) 

We note a few features of this method of implement- 
ing c-sign that will prove useful when using quantum error- 
correction to boost reliability. The first step in the implemen- 
tation of c-sign using \t' n ) is to perform BM n on the first qubit 
and the first n modes of \t' n ). With probability l/(n + 1), the 
measurement fails, and the first qubit is measured in the stan- 
dard basis with a known outcome. If this happen we do not 
attempt to complete the protocol so that the second qubit re- 
mains coherent. If the first measurement succeeds, the correc- 
tion step is applied. Note that as it only adjusts the phase, it 
commutes with the implemented operation. Next BM n is ap- 
plied to the second qubit and the third n modes of \t' n ). With 
probability l/(n + 1), this fails and the second qubit is mea- 
sured in the standard basis with known outcome. The outcome 
affects whether or not the first qubit experienced a sign flip, 
which can be corrected with a phase shifter to restore it to its 
original coherent state. 
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C. State preparation algorithm 

We have reduced the problem of implementing a near- 
deterministic c-sign to generating the state \t' n ). Clearly this 
can be done by first generating two copies of \t n ) and then 
applying the 0{n 2 ) c-sign operations. Since \t n ) is a uniform 
superposition of simple bosonic qubit states (the unary num- 
bers from to n), it is not hard to see how to construct a quan- 
tum network with 0(n) gates from an initial state accessible 
with single boson sources. Together this gives 0{n 2 ) gates, 
which can be implemented non-deterministic ally. Although 
n needs to be no larger th an so me constant (pessimistically 
no more than 25, see Sect. VIII , implying existence of tech- 
niques with asymptotically efficient resource usage), both the 
number of gates and the way in which non-deterministic gates 
are used needs to be improved. The remainder of this section 
is dedicated to this task. 

We begin by generating a variant of \t n ) with an ancilla 
bosonic qubit that contains parity information about the states 
in the superposition: 

n 

\tp n ) = £iim"-wii> n - j ' 

3=0 

x|(n - j) mod(2)) ai |(n - j + 1) mod(2)) a2 . 

(19) 

To prepare \t' n ) one can first prepare two copies of \tp n ), ap- 
ply c-sign to the two ancilla qubits, apply to each an- 
cilla qubit and measure them. There are four equiprobable 
outcomes to the measurement, which differ from the desired 
state \t' n ) only by signs, which can be undone by applying 
to the last n modes of one or both halves of the component 
states. This involves 0(n) individual steps, only one of which 
requires a non-deterministic element. 

To prepare state \tp n ), begin with the state (|l) n |0) n + 
v /n|l) n - 1 |0)|0)"- 1 |l))(|01) + |10)), which can be generated 
from a product of single boson states by applying beam split- 
ters to modes n and 2n and to the ancilla modes. Apply c-sign 
to modes 2n and 2n + 1. Let qi be the qubit encoded in modes 
I and n+l and a the ancilla qubit. The Tth step (0 < I < n— 1) 
consists of 

1 Conditionally on qubit being in the state 1 01) , ap- 
ply Bg t to qubit q n -i-i, with tan(^) = \fn — I — 1. 

2 Apply c-sign to modes 2n — I — 1 and 2n + 1. 

At the end undo the beam splitter operation on the ancilla 
qubit. 0(n) operations are required for implementing this 
algorithm and 0(n) of these are based on gates that we can 
only do non-deterministically. The steps of the algorithm can 
be implemented with c-sign, beam splitters and phase shifters 
using (for example) the methods of fll3|]. For example, the 
conditional beam splitter is obtained by conjugating a beam 
splitter on the target by c-sign and following that with the in- 
verse beam splitter. As this requires two c-sign operations, 
it may be more efficient to directly implement a conditional 
phase shift for phases other than —1 and then to conjugate 



this by a one qubit rotation to implement the required condi- 
tional beam splitter. If this is done by using an instance of our 
teleportation schemes (e.g. with n = 1), then the built in error 
detection may be exploited to achieve some error recovery. In 
particular, of the two failure modes for a gate implemented by 
teleportation, one results in a state that is essentially a \tpi) 
for I < k. This state can be used for improvements in the 
probability of success for the operations required when re- 
attempting the construction of \tp n ). 

The non-deterministic aspects of the method for preparing 
\tp n ) complicate the resource analysis. Let S(n) be the ex- 
pected number of elementary operations needed to prepare 
\tp n ). If we do not attempt to recover from failure, nor exploit 
improvements in the probability of success by using \tpk) for 

k < n, then the probability of success would be i° , lead- 
ing to S(n) = exponentially large. 

Using \tpk) recursively leads to a subexponential method 
as the following argument shows: Except for differences in 
rotation angles, the algorithm for preparing \tp n ) looks like 
the algorithm for preparing \tp n _\) followed by < C non- 
deterministic operations (for some constant C). Suppose we 
use 2C \tp ^f£) for these operations. Each of these states re- 
quires S(y/n) operations and their application has an addi- 
tional overhead of < D^/n. The success probability is is now 
1 - 1/y/n. Therefore, S(n) < (1 - 1/ \/n)-' 2C (S(n - 1) + 
2C{S{y/n) + Dy/n)). Under the assumption that S(ri) = 
fi(v^), S(n) < (1 + d/^)(S(n - 1) + C 2 S(v^)) for 
suitable choices of the constants. This implies that S(n) = 



D. Near-deterministic loss detection and non-destructive 
parity measurement 

Our constructions so far assume essentially perfect opti- 
cal gates and detectors. Unfortunately, loss of bosons is one 
of the primary error mechanisms. Since this implies "leak- 
age" errors for the bosonic qubits, scalability requires a non- 
destructive method for detecting when the modes supporting a 
qubit are no longer in a state associated with the bosonic qubit 
and returning them to a qubit state. Fortunately, our methods 
can be adapted to yield a non-destructive parity measurement 
for the space spanned by the states 1 00) , 1 1 0) , 1 1 ) , 1 1 1 ) . In or- 
der to implement such a measurement, it is sufficient to follow 
the procedure used for teleporting the c-sign gate, but using 
instead of \tp n ) the state 

IK) = E iOT> n -wiir- j - 

i,j:i+j=0(2) 

x|l} i |0} n - i |0) i |l) n - i . (20) 

(21) 

The measurement determines whether the number of bosons 
in the two input modes is even or odd, and the conditional state 
can be extracted from the output modes after suitable phase 
shifts. The state can be prepared using a simple variation of 
the method for preparing \t' n ), using only one ancilla qubit for 
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keeping track of the total parity. Measurement of the ancilla at 
the end yields either \p' n ) or the odd variant, which is equally 
useful. 

The above loss detection method is incomplete for detect- 
ing leakage in one respect: it does not detect when a mode 
has two or more bosons. However, it is sufficient for the ex- 
perimental procedures to guarantee return to the coding state 
with sufficiently high probability, that is, the loss need not be 
detected by the user. Although the information is lost in the 
process without that loss being detected, this is in principle 
sufficient for meeting the scalability requirements. The tele- 
portation methods described earlier for implementing various 
gates ensure that each output mode has at most one boson. By 
using them sufficiently frequently on both halves of a bosonic 
qubit, return to the state space where loss can be detected is 
thus adequately assured. 

Another application of the non-destructive parity measure- 
ment is to the teleportation of a bosonic qubit using the tra- 
ditional entangled state \E) = |0110) - |1001). The ini- 
tial state is (a |01) +otx\lQ))\E) on modes 1,2,3,4,5 and 
6, say. The Bell measurement requires measuring in the 
basis |0110) ± |1001) and |0101) ± |1010) This has been 
implemented experimentally with some probability of fail- 
ure [ p"4| , p"5| , Il6p. A non-destructive parity measurement on 
modes 2 and 3 determines which of the two pairs of basis 
states are present. The sign can then be determined by apply- 
ing two balanced beam splitters to modes 1, 2 and 3, 4 before 
measuring the four modes. 

To conclude this section, observe that the entanglement 
needed for teleportation can be produced locally, non- 
deterministically by performing the non-destructive parity 
measurement on modes 1, 2 of the non-traditional entangled 
state (|01X, 3 - |10\ ]3 )(|01^ 4 + |10^ 4 and accepting if the 
parity is odd. Of course, an appropriate entangled state is still 
used locally inside the parity measurement, but this one may 
be generated locally, perhaps using the techniques of this pa- 
per. The curious feature is that the distribution of the entan- 
glement requires only two independent bosons sent by means 
of beam splitters. The rest is taken care of by classical com- 
munication. 



VII. SINGLE BOSONS USING ACTIVE LINEAR OPTICS 

The need for non-linearities to create single boson states 
can be eliminated using weak squeezing and particle detec- 
tors [p7[]. A non-degenerate squeezer on modes 1 and 2 ap- 
plies the Hamiltonian s' 12 ) given in Eq. |j. When applied to 
the vacuum, the output consists of a superposition of states 
with identical numbers of bosons in the two modes. Under 
weak squeezing conditions, the relative amplitudes of states 
with a total of 2n bosons decrease exponentially with n. Thus 
a single boson can be produced in mode 1 conditional on de- 
tection of at least one in mode 2. The output state has arbi- 
trarily large overlap with the desired single boson state. The 
overlap is one in the limit of weak squeezing. If we can count 
bosons in mode 2, the reliability of the conditional output is 
even better. 



VIII. RELIABLE SCALABILITY 

That quantum computing is scalable in principle is a con- 
sequence of the accuracy threshold theorem jl8[ [T^ . |2C} Ejj ]. 
According to this theorem, the main requirement for efficient 
scalability of a physical implementation of quantum compu- 
tation is that the elementary operations can be implemented 
with a minimum accuracy, which is currently estimated to be 
less than 1 — 10 -4 [|22]]. Although achieving such accuracy 
may seem like a daunting task, in practice, errors behave much 
more predictably than assumed in the general analyses, and 
experience shows that at least some aspects of many quan- 
tum experiments are controllable with accuracies above the 
threshold. A typical example is the phase of RF pulses in nu- 
clear magnetic resonance. Such accuracies can be exploited 
to "boost" the accuracy of other gate parameters. 

Achieving the necessary accuracy for the methods intro- 
duced in this work is clearly not practical without additional 
work. For example, to implement a single two qubit gate with 
error rate < 10~ 4 would require 2 * 10 4 modes, half of which 
are initially loaded with single bosons. To get the necessary 
prepared state would require a very large number of tries in- 
deed. To avoid doing this we take advantage of the fact that 
ideally, it is always known when a gate fails. Thus, we begin 
by assuming that all optical operations (beam splitters, phase 
shifters, single boson state preparation, number state measure- 
ment) are perfect. Under this assumption, the threshold can be 
greatly improved by exploiting quantum erasure codes p3|], 
yielding an "erasure threshold". Given the special nature of 
the detected errors in our system, we are actually interested 
in an even more benign threshold Td- This results in quan- 
tum code enhanced systems, where qubits are highly protected 
in specially coded states. One can then estimate how errors 
in the operations propagate to errors in specific instances of 
such systems and apply the general threshold to obtain an es- 
timate of the minimum accuracy required. Most importantly, 
the overhead associated with using linear optics is directly re- 
lated to Td via the complexity of the necessary erasure code 
implementation. We pessimistically estimate Td to be well be- 
low .96, which implies that scalability is achievable using our 
teleportation protocols with size parameter n < 25. This still 
implies substantial overheads, albeit far from what might have 
been expected. The quantum code constructions that lead to 
our estimate of Td and more detailed resource analyses are 
in [H. 

IX. EXPERIMENTAL REQUIREMENTS 

A crucial question is to what extent the methods of the pre- 
vious sections can be implemented experimentally. Although 
large scale quantum computation is clearly still out of reach 
with current technology, many of the elements can be tested 
using existing equipment. For example, the non-deterministic 
c-sign requires only three modes and one ancilla single pho- 
ton. Similarly, the simplest methods for generating highly en- 
tangled states or for teleportation with parity measurements 
require reasonable overheads. 
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It is relatively easy to couple two optical modes at a single 
beam splitter. As subsequent beam splitters are added how- 
ever careful spatial and temporal mode matching is required, 
otherwise unwanted modes could be mixed in representing 
loss channels. Fortunately a four mode experiment could be 
achieved using only two spatial modes together with the two 
polarization degrees of freedom for each. The price would 
be to add a requirement for polarizing beam splitters, together 
with controllable polarizers (Faraday rotators), to enable flexi- 
bility in linear coupling and to separate the polarization modes 
for detection. 



The requirement of single photon sources is a more diffi- 
cult constraint to satisfy with available technology. To some 
extent this is mitigated by the conditional preparation of the 
non deterministic phase shift; while the time at which the gate 
is applied is random, we do know when it has occurred. One 
way around this difficulty is to use a pulse-trigger to indicate 
that temporally correlated photon pairs are incident to the de- 
vice. For example, in the case of the non-deterministic c-sign 
method, a three mode state with a total of three photons, to- 
gether with an ancilla mode in the vacuum state is required. 
The three mode state can be generated using a method similar 
to that used to generate a GHZ state [p5|] via a type II paramet- 
ric down converter with a suitable arrangement of polarizing 
beam splitters. To second order in the pump amplitude, the 
result is a four mode state with at most four photons in to- 
tal. Conditioning on a single count in a time window for one 
of the four outputs distributes three photons across the other 
three modes. Subsequent processing could then yield the re- 
quired input states to the non deterministic c-sign protocol. 
The considerable disadvantage of this approach is the very low 
probability of required states per run. Single photon sources 
however are not far away. The turnstile proposal of [ |26| l or the 
SAW method of [ 27 1 could yield pulses of light with a prede- 
fined maximum number of photons per pulse, synchronized to 
an electronic clock signal. Optical delay line methods could 
then be used to bring single photon pulses together at a par- 
ticular time and place. Of course single photon detectors with 
high quantum efficiency will be required in all cases. 



In an experiment it would be necessary to be able to distin- 
guish that the required controlled phase shift had been imple- 
mented in a successful run. The simplest way to do this is to 
use a double path interferometer with a controlled phase shift 
introduced in one arm only. The input signal state is then split 
at a 50/50 beam splitter at the input to the device, passes along 
two arms and is then recombined at an output beam splitter. 
In the absence of a phase shift the interferometer could be ad- 
justed to transmit the input state with unit probability at a sin- 
gle output port. When the phase shift is introduced however 
this probability would change as "which path" information is 
imprinted by the conditional phase shift. In essence the total 
experiment, including the conditional state preparation, is a 
kind of four photon coincidence experiment. 



X. DISCUSSION 

The ability to implement quantum computation with linear 
optics and particle detectors realizes the dream of computing 
with non-interacting particles. The only particle interactions 
occur implicitly in the detectors and result in particle destruc- 
tion. The present work shows how to do this for bosons, but 
similar techniques work for fermions, which have the prop- 
erty that the creation operators satisfy anti-commutation laws. 
This greatly improves on previous methods for implementing 
quantum networks with linear optics [^, BQ], which re- 
quire an exponential number of modes to represent the state 
space of n qubits. It solves a problem we first learned from 
Paul Kwiat, who asked what the computational power of pre- 
pared entanglement in optics is. Other studies of quantum 
computation with harmonic oscillators or continuous vari- 
ables have depended on non-linear effects either for quantum 
gates |[| [H]. [32|, [33|] or for state preparation [Q|, 

LOQC is related to the observation that in the usual model 
of quantum computation, it suffices to be able to implement 
operations in the normalizer of the Pauli group, provided 
one non-stabilizer state can be generated, such as the state 
cos(7r/8)|0) + sin(7r/8)|l) [^0|, or a non-stabilizer measure- 
ment can be performed, such as that of the Hadamard opera- 
tor H. Linear optical elements generate the normalizer of the 
continuous version of the Pauli group. This group is given by 
products of operators of the form e -4 '( a+at ) and e -t( - a ~ at - ) , 
which correspond to translations in phase space. From this 
point of view, the vacuum state preparation and particle de- 
tectors serve the purpose of the non-stabilizer state and mea- 
surement, respectively. 

Even though LOQC is no more practical given current tech- 
nology than any of the other proposals, many of the basic el- 
ements can be tested, and there are numerous other applica- 
tions in quantum information processing. For example, our 
methods can be used to generate complex entangled states 
with higher efficiency than is possible by using sequential 
down-conversion sources. The non-destructive parity mea- 
surement can be used for implementing standard quantum 
teleportation essentially unconditionally. In general our meth- 
ods yield complex quantum non-demolition measurements. 
Since quantum communication is substantially less demand- 
ing than quantum computation (as demonstrated by its much 
better thresholds [[34]]), our methods are closer to being prac- 
tical for implementations of quantum channels, particularly 
when relatively low success probabilities are acceptable, such 
as in quantum cryptography. Note that quantum optics is 
likely to dominate the important communication sector of 
quantum information processing regardless of how quantum 
computation is finally implemented. 

Another approach to implementing quantum computation 
optically is to use non-linear optical elements. Assuming that 
such elements are available, though perhaps with high losses, 
it may be possible to combine techniques from LOQC with 
these elements to significantly reduce the introduced errors. 

More work is required to optimize the resources required 
for LOQC. Clearly the suggestions made here are just the be- 
ginning and we believe that significant improvements are pos- 
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sible. From a theoretical point of view it would be nice to 
establish lower bounds on the complexity of various elemen- 
tary gates depending on the desired probability of failure. 
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